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Abstract 

We present an analysis of the m 2 -corrections to Cabibbo-suppressed r lepton de- 
cays employing contour improved resummation within an effective scheme which is 
an essential new feature as compared to previous analyses. The whole perturbative 
QCD dynamics of the r-system is described by the /3-function of the effective cou- 
pling constant and by two 7-functions for the effective mass parameters of the strange 
quark in different spin channels. We analyze the stability of our results with regard 
to high-order terms in the perturbative expansion of the renormalization group func- 
tions. A numerical value for the strange quark mass in the MS scheme is extracted 
m s (M T ) = 130 ± 27 exp ± 9 t h MeV. After running to the scale 1 GeV this translates 
into m s (l GeV) = 176 ± 37 cxp ± 13 th MeV. 

PACS: ll.10.Hi, 12.38.-t, 13.35.DX, 14.65.Bt 
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1 Introduction 



The r-system offers the possibility to confront QCD with experiment in the low en- 
ergy region. The high precision of experimental data and good accuracy of theoretical 
results make r-physics an important testing ground for QCD j[J ^ ||. Theoretically 
the observables of the r-system are related to the moments of the spectral density 
of a correlator of hadronic currents which can be reliably calculated within pertur- 
bation theory [|], [5], ||. Therefore, the r-system observables were extensively studied 
during the past few years within the operator product expansion (OPE) which is a 
general approach to analyzing the properties of current correlators 0, [| [TD|, |TT 



The perturbation theory (PT) series in QCD appear to be asymptotic and the ul- 
timate accuracy they can provide depends on the concrete numerical value of the 
expansion parameter - the strong coupling constant a s (E) at a relevant energy E. 
This limits the theoretical accuracy which can be obtained within the finite order per- 
turbation theory (FOPT) analysis. In the case of the r-system the strong coupling 
constant a s (M T ) is not small at the scale of the r lepton mass M T which can lead to 
an asymptotic growth of terms of the perturbation theory series already at a rather 
low (third-fourth) order of PT expansion. Judging from the analysis of the moments 
of the hadronic spectral density in the finite energy interval (0, ) within finite 
order perturbation theory there are strong indications that the ultimate theoretical 
accuracy for the r lepton decay observables has already been reached at next-to-next- 
to-next-to-leading order (N 3 LO) or a z s which is the highest order of PT expansions 



presently available |T2|- The convergence behavior of the perturbation series for the 
r lepton observables depends on the region of the spectral density which is being 
probed: if the low energy region is suppressed (as for high moments of the spectral 
density) the asymptotic limit of the series moves to higher order terms. The expan- 
sion of the correlator in m? s makes the explicit convergence of the PT series for the 
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coefficient functions of consecutive m 2 s corrections slower JTjJ. The reason for such 
a behavior is quite obvious - higher order terms of m 2 s expansion of the correlator 
are more sensitive to the low energy region of integration in Feynman diagrams that 
is not described by perturbation theory. Thus, the PT expansions for observables 
in the r-system seem to be at the edge of asymptotic growth in the N 3 LO. At the 
same time, the present accuracy of experimental data for some observables of the 
r-system is already comparable with the ultimate theoretical accuracy reachable in 
FOPT 0, [3] . This raises the problem of obtaining more precise theoretical formulas. 
Higher order terms of FOPT (though they are very desirable and provide additional 
information) will not give more precise results for PT series if the asymptotic limit 
is already reached and must simply be discarded in FOPT applications. To catch 
up with improving experimental accuracy it is necessary to find a way of interpreta- 
tion of the perturbation theory series that allows one to reach a theoretical precision 
comparable to that of the experimental data. One possibility to extract numerical 
results from a perturbation theory series which is explicitly divergent at finite or- 
ders is to apply a resummation procedure (see e.g. This is more sophisticated 
than just summing the consecutive terms of the perturbative expansion up to some 
finite order but it requires some knowledge (or assumptions) about the behavior of 
infinite number of terms of the expansion. The choice of a resummation procedure 
is not unique and there are many ways to resum or improve the convergence of an 
asymptotic series e.g. J15, [U|. We think that there are two important criteria for 
the choice of an appropriate resummation procedure for perturbation theory series of 
physical observables: the renormalization group structure of the PT series should be 
respected and the definition of the effective parameters used for the description of the 
observables should be physically motivated. Led by these two criteria we use contour 
improved perturbation theory (CIPT) (10], |T7| (see also [[H|) in an effective scheme 
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(not MS) to resum perturbation theory contributions in all orders. Within CIPT, the 
procedure of resummation is, however, not unique. A freedom in the choice of dif- 
ferent renormalization schemes still remains within CIPT, which will affect the final 
numerical results |1| . We choose an effective scheme which we consider to be natural 
and the simplest one for the r-system. Note that the choice of an adequate scheme 
is also dictated by the way the system is to be described. In the FOPT analysis of 
the moments, the contour integration is completely (mathematically) equivalent to 
the integration along the physical cut. Therefore, the basic object for FOPT analysis 
is the PT spectral density of the correlator which determines the effective coupling 
in this procedure |f[2|| . In CIPT, the correlator is the basic perturbation theory ob- 
ject which is naturally normalized in Euclidean domain. Within the effective scheme 
approach all PT corrections to the correlator are absorbed into the definition of the 



effective parameters of the system |2(J, glj, [22|, |24], P5[ . In the massless case the 
effective coupling constant is the only relevant parameter. If the m s -corrections to 
the correlator are included then two additional parameters m 2 and m 2 related to 
mass corrections to different spin structures of the correlator should be considered. 
These parameters are perturbatively related to a finite strange quark mass. In the 
massless case the effective scheme resummation analysis along these lines was done 



in [fL9|| . In this paper the resummation analysis in the effective scheme is extended to 
m 2 corrections. An analysis of the m 2 corrections within CIPT in the MS scheme was 
previously performed in refs. |26|, p7]. Our results of resummation within the effective 
scheme approach confirm the main conclusions of the paper [26[ and are consistent 
with the results of ref. [[27]] though the direct comparison is a bit complicated because 
we use a different parameterization of the relevant observables. We also use differ- 
ent moments for the final determination of a numerical value for m s as compared to 
ref. p?| - the detailed discussion of our choice is given in the text. The uncertainty 



4 



of our new results is smaller than the conservative error bars given in ref. on 
the base of an analysis of the explicit convergence of PT series. In order to get an 
understanding of the reliability of our present procedure we evaluate the stability 
of our results with regards to higher order corrections to the renormalization group 
functions. 

The paper is organized as follows. In Sect. 2 we set up the stage for resum- 
mation in an effective scheme. We define our basic quantities and calculate their 
respective anomalous dimensions (beta-functions) that determine the running along 
the integration contour. The known results of PT calculation in the MS scheme are 
used but the expressions for the corresponding renormalization group (RG) functions 
of effective parameters are new and specific for the effective scheme defined in the 
paper. In Sect. 3 we present the results of theoretical calculations of the moments 
of the spectral density of r decays in the effective scheme. The PT expressions are 
given in the effective scheme which is a new result. We introduce the effective mass 
parameters for the channels with spin 1 and spin and give explicit formulas ex- 
pressing these parameters through the standard MS mass m s in a resummed form. 
Some nonPT corrections are taken at the leading order only and are therefore scheme 
independent: they coincide with the results already published in the literature (see, 
e.g. ||27|| ). In Sect. 4 we extract a numerical value for the strange quark mass from 
experimental data and compare our results with conclusions of earlier analyses. Sect. 
5 gives our conclusion. In an Appendix the RG functions of the effective scheme are 
explicitly given through the coefficient of the MS-scheme RG functions known from 
the literature. 
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2 Resummation in an effective scheme 



In the first step of our analysis of the r-system we define an effective scheme in which 
all higher order PT corrections to the correlator of hadronic currents are absorbed into 
the effective coupling and two effective mass scales, or two coefficient functions of mass 
corrections. If such a scheme is used for describing the r-system the only perturbative 
objects are one effective /3-function for the coupling constant and two effective 7- 
functions for the coefficient functions of mass corrections. Using these three functions 
we determine the evolution of the effective coupling and of the two coefficient functions 
of mass corrections on the contour in the complex plane of squared momentum. 
Within our procedure the (5- and 7-functions are treated as exact functions which 
is a standard understanding of renormalization group (RG) summation. Having the 
explicit solutions for the running coupling and the mass coefficient functions in hand 
we determine observables such as the moments of the decay rate simply by integrating 
the coefficient functions with a weight function specific to a chosen observable. 

2.1 Definition of an effective scheme 

The basic theoretical quantity for describing the r semileptonic decays is the correlator 

of two hadronic currents 

r N 
IW?) = i J dxe^(T 3ll (x)jt(0)) = ^(g^IT g (g 2 ) + g^U g (q 2 )) (1) 

with jn(x) = M7 M (1 — 75)5- Here Yl q {q 2 ) and Ti g (q 2 ) are scalar invariant functions, 
N c — 3 is the number of colors in QCD. We work within QCD with three light quarks. 
The correlator is normalized to unity in the leading (parton model) approximation 
with massless quarks. By expanding Ti q (q 2 ) and Ti g (q 2 ) in terms of powers of m 2 /q 2 

and keeping only the first order term in this expansion one has 

2 

in 

U q (q 2 ) = U(q 2 )+3^U mq (q 2 ), (2) 
6 



II 9 (g 2 ) = -g 2 n(g 2 ) + \m 2 s Il mg { q 2 ) (3) 

where IT(g 2 ) is an invariant function already known from the mass zero case. The 
functions Tli(Q 2 )(i = m q ,m g ) with Q 2 = —q 2 are computable in perturbation theory 
in the deep Euclidean region Q 2 — > oo. The results of perturbation theory calculations 
for the correlator given in eqs. (^J|) were obtained in refs. ^9], 30, 31] and have 
already been used in the FOPT analysis [E 



We define new effective quantities a, m 2 , m 2 such that all information from per- 
turbation theory calculations for the functions U(Q 2 ), H m q,mg(Q 2 ) is absorbed into 
the evolution of these new quantities which is determined by the effective (3- and 
7-functions. For the mass corrections we introduce two different coefficient functions 
(we sometimes call them the effective mass parameters) because the correlator in 
eq. (HI) is not transverse if corrections of the order m 2 are taken into account. Our 
definitions of effective quantities a, m 2 m 2 are 



Q 2 ^(Q 2 ) = l + a(Q 2 ), 



d 

dQ 2 
d_ 
dQ 

,2/ii^-2\tt ir& 



-<(M^)Q 2 ^-U mg (Q 2 ) = m 2 g (M 2 )C g (Q 2 ) , (4) 



m 2 s (M 2 )Tl mq (Q 2 ) = m 2 q (M 2 )C q (Q 2 ) . 

Here C q , g (Q 2 ) are coefficient functions of mass corrections. They are conveniently 
normalized by the requirement C q ^ g (M 2 ) = 1. 

In terms of the MS scheme quantities a s = a s (M 2 ) and m s = m s (M 2 ) the effective 
parameters in eq. @ read 

a(M 2 ) = ^ + k 1 (^)\k 2 (^Y + k 3 (^-Y + 0(a 5 s ), (5) 

7T V 7T / V 7T / V 7T / 

m 2 (M 2 ) = ^(^(i + ^ + ^^ + ^^' + ^ai)), (6) 
m 2 q (M 2 ) = m 2 (M 2 )(l + ^ + ^(^) 2 + A; g2 (^) 3 + 0(^)). (7) 

Numerical values for the coefficients k%, k q 2 are unknown though their estimates within 



various approaches can be found in the literature. We have written explicitly first 
coefficients for mass corrections, k g o = 5/3, k q o = 7/3. Further references concerning 
numerical values for the known coefficients ki^, kgi,g2, k q i are given in Appendix. 



2.2 Running of the effective coupling a(Q 2 ) and 
the mass coefficient functions C g g (Q 2 ) 

The behavior of the effective coupling a(Q 2 ) and the coefficient functions of the mass 
parameters C qt9 (Q 2 ) is determined by the effective beta and gamma functions. The 
defining RG equations for the evolution of these effective quantities are 

Q 2 ^ a (Q 2 ) = (8) 



Q 2 ^-2° 9 (Q 2 ) = 2 lg (a)C g (Q 2 ), (9) 



Q2 W 2 ° qm = 2 ^ C ^- ( 10 ) 
The RG functions /3(a) and 7 g , g (a) describing the evolution of the parameters in the 
effective scheme can be expressed through the MS scheme RG functions (the standard 
/5-function of the coupling constant and the mass anomalous dimension 7) using the 
relations in eqs. (f|,[5],|6],|7|) . The explicit formulas of such a RG transformation are given 
in Appendix. Up to the relative order a 3 , the RG functions f3(a s ) and 7(a s ) in the 



MS scheme have been calculated in refs. |32|, B3|, |34|. Under the RG transformation 



eqs. (§,0,0) the MS scheme j3- and 7- functions are transformed into the effective (3- and 
7- functions (§,|||10]). Note that first two coefficients of the /3-function /3 , /?i and first 
coefficient of the 7-function 70 are invariant under RG transformations. Numerically 
the RG functions of the effective parameters describing the r-system up to the order 
m 2 s are given by 

(3{a) = -a 2 (2.25 + 4a + 11.79a 2 + a 3 (-76.36 + 4.5fc 3 )) , (H) 



7 P (a) = -a (l + 4.027a + 17.45a 2 + a 3 (249.59 - fc 3 )) , (12) 
7,(a) = -a (l + 4.78a + 32.99a 2 + a 3 (-252.47 - k 3 + 3.38fc, 2 )) • (13) 

For the effective coupling in eq. (|J) we use the numerical value a(M 2 ) = 0.1445 



35 1 extracted from the r decay rate into nonstrange particles within the effective 
scheme resummation procedure described in ref. This value corresponds to the 
MS scheme value of the coupling constant a s (M 2 ) = 0.343 that is a bit larger than 
the most recent result obtained in FOPT analysis of the r decay rate into nonstrange 
particles |3"5" |. The PT series for the effective /3-function in eq. fllTD explicitly converges 



well at the numerical value a(M 2 ) = 0.1445. If the coefficient k 3 lies in the range 
< A?3 < 50 (which is a conservative estimate based on a number of predictions) the 
a 3 coefficient in eq. (|11]) will not be extremely large but nevertheless the /3-function 
(|TT|) shows asymptotic growth in the N 3 LO for k% > 35. The 7 9 -function in eq. (12 



behaves worse than the /3-function but still the explicit convergence persists up to 
the NNLO at a(M 2 ) = 0.1445. The N 3 LO correction will show asymptotic growth 



for values of k$ smaller than 129. The 7 g -function in eq. (|T3|) has already shown an 
asymptotic growth in the NNLO which will limit the precision of our results. The 
(3- and 7 3 -function will show asymptotic growth in the N 3 LO because no choice of k 3 



can make them convergent simultaneously. This confirms the conclusions of ref. |p~2| 
where the asymptotic growth for the FOPT expressions of the moments of the spectral 
density has been found in N 3 LO independently of the choice of the numerical value 
for k 3 . In other words, the evolution of the effective parameters for the r-system is 
too different to be handled by the FOPT expressions for the RG functions in the 
N 3 LO. 

The running of the effective coupling and the mass coefficient functions C qt9 (Q 2 ) 
along the contour in the complex plane of momentum squared is determined by the 
renormalization group equations. It is convenient to choose a circular contour in the 
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complex Q 2 -plane and to parameterize it by the relation Q 2 = M 2 e 1 ^, —it < <p < n 
which leads to the differential equations for the coupling constant 

-i4w = ^)), a(0 = 0) = a(M T 2 ) (14) 

and for the coefficient functions C n {(j>) 

- i^C n {<f>) = 2 7n (a(0))C n (0) , C n (0 = 0) = 1 (15) 

with n — q,g. The solution to the differential equation for the running mass in terms 
of the coefficient functions C n {4>) in eq. fllSD can be expressed through the integral 

C n {<j>) = exp hi j\ n (a( X ))dxj . (16) 

The initial values for a(4>) and C n (4>) are fixed at Q 2 = M 2 or cf) = 0. All corrections 
stemming from higher order perturbative terms are absorbed into the (3- and 7- 
function coefficients if the effective scheme is used as it is defined in eq. (f|). The 
solutions to the differential equations for the coupling a(0) (eq. (|T4|)) and for the 
coefficient function C g (0) (eq. (|T5|)) are shown in Figs. (]1]J2|)- The effective coupling 
a(4>) does not change much uniformly on the contour when higher order corrections 
of the /5-function are included. Especially the change from NLO to NNLO is rather 
small. The behavior of the coefficient function C g {4>) is rather similar to that of the 
effective coupling a(4>) as one can expect from the structure of PT series for j g (a) in 
eq. ( flit) in comparison with the /^-function in eq. (0). The function C g (<j)) appears 
to converge uniformly in the interval \(f)\ < it when going from NLO to NNLO. The 
function C q (<f>) (Fig. (0)) does not converge uniformly on the contour because of the 
slow convergence of the 7 g -function (eq. (0)) in consecutive orders of PT expansion. 



10 



-3-2-10123 -3-2-10123 

tp rad tp rad 



Figure 1: Running of the effective coupling a(0) on a circular contour in the complex 
plane calculated at LO, NLO and NNLO (left: real part; right: imaginary part) 




-3-2-10123 -3-2-10123 



tp rad tp rad 

Figure 2: Running of the mass coefficient function C q (<f)) on a circular contour in the 
complex plane calculated at LO, NLO and NNLO (left: real part; right: imaginary 
part) 

2.3 Resummation 

We use the direct generalization of the resummation procedure described for the 
massless case in ref. [[19[]. We treat the renormalization group functions f3(a) and 7 n (a) 
in any fixed order of PT (LO, NLO and NNLO) as exact functions and solve exactly 
(though numerically) the differential equations (0,0) for the running parameters 
a(0) and C ? 5 (0) on a circular contour in the complex plane. 
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Having the solutions for the running parameters we calculate numerical values 
for the theoretical expressions of the observables in the r-system. The observables 
of interest are expressed through the moments of the spectral density of the basic 
correlator given in eq. ([!]). For the massless part of the correlator H(q 2 ) given in 
eqs. we define moments of the spectral density by the relation 

mm = (B + 1) ^/n (9 ")(^)"g 

= l + (-l)"-!- f e«**'l*a(^+i f O (0)#. (17) 

For the mass correction related to the g-part of the correlator the expressions for the 
moments are very similar. One has 

M 9 (n)= [ —^ e «°+i)*c a (<t>)d4>+ 7r C s (0)#. (18) 
For the mass correction in the g-part of the correlator we use the definition of physical 



moments given in 12 . The corresponding expression on the contour reads 

Mf (n) = (-l) n+1 ^ f C ? (0)e^#. (19) 

The moments of the massless part M(n) and the mass corrections M q>g (n) are the 
basic objects which can be calculated theoretically. Our qualitative conclusions about 
the stability of the running parameters due to the higher order corrections to the f3- 
and 7-functions are confirmed by the behavior of the moments. The moments of the 
massless spectral density given in eq. fllTf) are stable with regard to inclusion of higher 
order corrections of the /^-function. This is a direct consequence of the /3-function 
behavior in consecutive orders of PT along the integration contour plotted in Fig. |l[ 



The moments of the g-part given in eq. ( 18|) are less stable with regard to inclu- 



sion of higher order corrections of the 7 9 -function than the moments of the massless 
part. The moments of the g-part in eq. ([19]) behave worse in higher orders than the 
moments of the g-part which reflects the slower convergence of the 7 g -function in 
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eq. fll3"|) as compared to 7 9 -function in eq. fll2|) . In the massless and in the g-part 
the values of the moments are dominated by the second integral in eqs. ([T7|JT8|) be- 
cause the values of the integrals with oscillating integrands in eqs. (|T7| , p!8|) decrease for 
high order moments. The coefficient functions C q ^ g (Q 2 ) depend on the 7,^-functions 
exponentially according to eq. (|16|) and they are very sensitive to the convergence 
of the 7g i9 -functions. The exponential dependence of the mass coefficient functions 
on RG functions is an essentially new feature of the analysis of the m? s corrections 
as compared to that of the massless part of the r lepton decay rate. Such strong 
dependence requires more accurate treatment of PT series for RG functions. 



3 Determination of m s from r decays 

One of the important aims of the analysis of Cabibbo-suppressed r decays is the 
extraction of a numerical value for the mass parameter m s . Various observables 
can be used for this purpose. Here we consider the r lepton total decay rate into 
hadrons. Experimental data for Cabibbo-suppressed hadronic r decays are not yet 
very precise. Contrary to the experimental data the theoretical expressions for the 
relevant observables related to the moments of the spectral density of the correlator of 
the hadronic currents are known rather precisely in a sense that the corresponding PT 
series within OPE are calculated in high (next-to-next-to- leading) orders of expansion 
in the strong coupling constant. However, this theoretical accuracy is quite formal, 
because the explicit convergence of PT expressions for the moments is rather slow 
and the PT series has to be resummed in order to obtain numerical results. 

The theoretical expression for the m 2 s corrections to moments (k, I) of the differ- 
ential decay rate is given by the contour integral in the complex g 2 -plane 
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j^A kl + ^B kl j = -6^ (u q A kl + u g B kl ) = -6^F kl (20) 

where the superscript (k, I) denotes an integration with additional weight factors 
which suppress the high (k > 0, 1 = 0) and low (k = 0,1 > 0) energy region. The 
finite order PT expressions for the coefficients F k i can be found in refs. [0, |26j . The 



analyses performed in ref. [ 12 , pB| in FOPT and CIPT within resummation in the 
MS scheme have demonstrated that the explicit convergence of the PT series for the 
moments is slow. A close conclusion about the convergence of PT series in the MS 



scheme was given in ref. |27|. Therefore, in the present paper we use the contour 
resummation procedure in the effective scheme and calculate the relevant coefficients 
F k i in a closed form. 

We use the (0, 0) moment as our best theoretical estimate for comparison with 
experiment. Other moments are briefly discussed but are not used in numerical anal- 
ysis for the strange quark mass determination. The choice of optimal moments is 
discussed in detail later. In brief, the reasoning is based on analyzing which energy 
regions saturate the respective moments. Indeed, to choose a special linear combina- 
tion of moments is equivalent to integrating the running parameters with a special 
weight function. Whether the corrections to a given linear combination are large or 
small can depend strongly on the integration region in which the integral of the spec- 
tral density is saturated. The change of the running parameters with higher order 
corrections is big in the region close to the physical cut (at = ±7r) and small in the 
deep Euclidean region where the functions are fixed by their starting values. If some 
observables have the regions <fi = ±7r strongly suppressed they have no contribution 
from the perturbative running and, therefore, are very stable with regard to the inclu- 
sion of higher order terms of RG functions. The moments of the type (k, 0) represent 
such observables for large k. However, these moments have the contributions from 
large energy region strongly suppressed and receive relatively large contributions from 
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the low energy region that makes them rather nonperturbative. Therefore, the use 
of the moments (k, 0) with large values of k is not under reliable quantitative control 
within OPE and perturbation theory, though their use is favorable from the point of 
view of precision of experimental data. We elaborate on this point later. 



3.1 Relation of the MS mass to effective mass parameters 

Within the effective scheme approach the r-system is described in its own terms with 
the mass parameters m qyg . In order to obtain a numerical value for m s which can be 
compared with other determinations we express the natural mass parameters of the 
r-system vn q g through the MS mass parameter m s . We emphasize that this is only 
done for purposes of comparison. In principle, observables in the r-system are best 
described by their internal parameters a, m q and m g . Relations between observables 
within the r-system can be found without any reference to the standard MS scheme 
parameters. 

In general, the perturbation theory expression for a given observable in a given 
renormalization scheme is parametrized by a mass scale and by coefficients of the 
evolution (RG) functions. In the massless case these are the scale parameter A and 



the perturbation theory coefficients of the /3-function (e.g. |36|, 37j). In the massive 



case there are in addition the invariant mass M to be defined in eq. ([H]) and the 
coefficients of the 7-function describing the evolution of the running mass. As in 
the case of the scale parameter A, the invariant mass M can be defined in different 
ways. The concrete definition may be fixed by a given asymptotic behavior at large 
momenta. This is the way the standard scale parameter A^jg is fixed. We define the 
invariant mass Mina/i independent way by writing 

Note that M is renormalization group invariant. That means that if m(/^ 2 ) is redefined 
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by some RG transformation the change is absorbed by the corresponding change of 
the 7-function so that M remains invariant up to the order in the coupling which 
has been taken into account. The RG invariance of M can be used to relate running 
masses defined in different renormalization schemes. After squaring eq. we find 
the relation between two mass definitions in different schemes with the 7-functions 
7(a) and 7'(a) expressed in terms of the same coupling a. In our analysis we use the 
effective coupling a as it is defined in eq. @ and relate m\ (n = q, g) to m' 2 by 

«W ) = m a {p )expj -2^ — d£j . (22) 

Here m' 2 is the standard MS scheme mass but with the evolution function 7' (a) 
expressed through the effective coupling. Eq. (|22|) relates two mass parameters de- 
fined in different renormalization schemes but expressed through the same coupling. 
Therefore the standard 7-function of the MS mass should be rewritten in terms of 



the effective parameter a before the use in eq. (|22|). It is also possible to use eq. ( pT|) 
written in the MS scheme and in the effective scheme. Then one determines the 
relation between m 9)9 and m s by eliminating the invariant mass M directly 



In our case the effective couplings a n , n = q,g are equal to the same coupling a 
for both masses m q>g . The two procedures described above lead to close numerical 
values for the coefficients ui q%g relating the effective mass parameters m q ^ g to the MS 
scheme mass m s . The difference of numerical values for u q and uj 9 obtained from the 
two procedures described by eqs. (|22|) and ( p3[) turns out to be less than 3% which 



is the residual scheme dependence of the results. Eqs. (|22"| , [23D allow us to express 
the internal parameters m q , g through the standard MS scheme parameter m s (see 
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eqs. and (j9|, PH) ). Finally we find for the coefficients u q>g relating the effective 

mass parameters m 9 5 to the reference MS mass m s 



m q = uj q m s , m g = u g m s , 



Ug = 1.73 ± 0.04, oj g = 1.42 ± 0.03 . (24) 

The numerical values for the coefficients u q>g are not close to unity which shows that 
perturbation theory corrections for observables in the r-system are rather large in 
the MS scheme. The FOPT expansion for the coefficients u qtg has been obtained 
in ref. [HJ. This expansion converges slowly which forces us to use the exact RG 



conversion given in eqs. (E2LE3p 



3.2 Power corrections from dimension D = 4 
condensate terms 

For the determination of m s one needs not all D — 4 condensate corrections to 
the theoretical expression for the r lepton decay rate but only those that enter the 
difference 

pfcZ TDkl 
\ v ua\ | »«s I 

Here Rf s=Q x is defined as 

and dR TS= o,i/ ds is the differential r-decay rate into the final hadronic states with the 
strangeness 0, 1 and the energy y/s. In the theoretical expression for the difference 
we neglect terms of the order mf/M^, set the u- and ci-quark masses to zero, and 
retain only the most important term linear in m s (cf. refs. [^, |2"7|]). Within OPE the 
coefficient of this term is given by the quark condensate. The final result reads 

ffl" = N cSew (6^F t , - 4^||n,) (27) 
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where N c = 3, m s = m s (M^), and = 1-0194 describes electroweak corrections 



39 . The coefficient function of the D = 4 local operator m s ss is taken in the 
leading order of perturbation theory expansion within OPE. In this approximation 
for the coefficient function the quantities multiplying the quark condensate are 
given by the expression 

T kl = 2 {8 lfi {k + 2) - 8 l}1 ) . (28) 
The numerical values for the first few coefficients T k i read 

T o = 4, T 10 = 6, T 20 = 8, T 01 = -2, T n = -2 . (29) 

These results agree with the leading order expressions for the coefficients given in 
ref. f27|. We have set up all ingredients necessary for the evaluation of eq. (|27|) . 
Table [l] gives the coefficients of the q- and g-mass parameters in eq. (^) . 



4 Numerical analysis and the choice of moments 

Having theoretical expressions for all moments we have still to optimize the choice 
for comparison with experiment and extraction of the numerical value for the strange 
quark mass. As was mentioned above the theoretical expressions for the moments 
(0,/), I > are more reliable from the point of view of perturbation theory than 
those for the moments (k, I) with nonzero k (the detailed discussion of this point is 
given below). However, experimental precision is worse for the moments (0,/) with 
large / because such moments are saturated by the contributions of many-particles 
hadronic states which is difficult to measure (see Table 0) . Note that some many- 
particles hadronic state contributions in the experimental data (for instance, K4it 
contribution) are represented by a result of Monte-Carlo simulation rather than di- 
rect measurements. Therefore, we use the moment (0, 0) as our best choice from both 
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(k,l) 


ALO 
A kl 


4NLO 
A kl 


/1NNLO 
kl 






£>NNLO 


(0,0) 


1.361 


1.445 


1.434 


0.523 


0.601 


0.625 


(1,0) 


1.568 


1.843 


1.976 


0.441 


0.552 


0.601 


(2,0) 


1.762 


2.282 


2.646 


0.390 


0.530 


0.607 


(0,1) 


-0.207 


-0.398 


-0.542 


0.082 


0.050 


0.025 



Table 1: Coefficients of eq. (20) 



experimental and theoretical point of view. The theoretical expression for this mo- 
ment exhibits a rather good convergence in consecutive orders of perturbation theory 
and the accuracy of experimental data for it is still acceptable in comparison with 
higher (0, I) moments. Note that the perturbative convergence is the main concern 
of the theoretical analysis in both massless and massive cases. In the massless case 
the nonperturbative corrections are small if factorization is used for the four-quark 



condensates pj BO, 41 . 



4.1 Numerical value for the strange quark mass 

The coefficients Am are related to the g-part of the correlator. This part contains 
contributions from spin and spin 1 particles. The spin piece is prone to possible 
nonperturbative contributions of (direct) instantons and perturbation theory expan- 
sions are expected to break down in low orders in this channel. However, this is only a 
general expectation without strict quantitative estimates of applicability of PT. The 
coefficients are related to the g-part of the correlator. The g-part contains only 
contributions of spin 1 particles. Nonperturbative contributions of (direct) instantons 
are forbidden in this channel by symmetry considerations. The coefficient Fki of m 2 s 
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term in eq. ( p7|) is combined according to eq. ( |20| ) 

F kl = u q A kl + u g B H . (30) 
As for the linear in m s term in eq. ( p7|) we calculate its numerical coefficient using a 



phenomenological value for the quark condensate. We use the relation [42, 43, 44 



(ss) = (0.8 ±0.2) (uu) (31) 

and the numerical value (uu) = —(0.23 GeV) 3 which coincides with the standard 
value (see e.g. f2~7fl ). Substituting all necessary quantities into eq. (|27j) we arrive at 
the defining equation for X = m s j (130 MeV) 



2 

-,kl 



;vsb (ssv)' 6R " = x m ' ■ x + ° 936 ' Tu) ■ (32) 

The dimension-four term contributes appreciably to the total theoretical result for 
different moments. For the moments (k, 0) its relative contribution increases with k 
for the first few moments. We use only the moment (0, 0) for which the dimension- 
four contribution gives about 16% of the total result. The coefficient function of 
the dimension-four contribution converges well in the perturbative expansion within 
OPE. We do not take the PT corrections to the coefficient function of the dimension- 
four contribution into account (see e.g. pTfl ) because of the large uncertainty in the 
numerical value of the strange quark condensate in eq. (|31~1) . 

For extraction of the numerical value for the strange quark mass we use the ex- 
perimental data obtained by the ALEPH collaboration 0]. The results for m s (M^) 
from different moments are given in Table |2|. For the determination of m s we use 
only the moment (0, 0) as the most reliable one from the perturbation theory point 
of view. A detailed discussion of the justification for our choice is given later. 

The final relation for determining the central value of the strange quark mass from 

the data on Cabibbo-suppressed r decays reads 
1 f M \ 2 



N C S EW V130 MeV 



J (ffiu^exp = 24.1 = x (20.2X + 3.74) . (33) 
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(M) 




m s (M r 2 ) MeV 


(0,0) 


0.394 ±0.137 


130 


(1,0) 


0.383 ±0.078 


111 


(2,0) 


0.373 ±0.054 


95 



Table 2: Results for m s (M T ) obtained from different moments of 5R T 

The result is X = m s /(130 MeV) = 1.00 . . . with the accuracy of two decimal places. 
This leads to our final prediction for the strange quark mass at M T 

m s {M 2 T ) = 130 ± 27 cxp ± 3 <s - s> ± 6 th MeV . (34) 

Note that this result is obtained with the numerical value for the effective coupling 
a(M^) = 0.1445 extracted from the r decay rate into nonstrange particles within 
the effective scheme resummation procedure described in ref. The reference 

numerical value of the MS scheme coupling constant is a s (M 2 ) = 0.343. 

We also give the value for m s (l GeV) obtained from m s (M 2 ) in eq. fl3"4| ) after 
four-loop running in the MS scheme 

m s (l GeV) = 176 ± 37 cxp ± 4< 5s> ± 9 th MeV . (35) 

The numerical value of the strange quark mass at the scale 1 GeV depends on the 
way the evolution is performed because of the truncation of the PT series for the RG 
functions. This difference reflects the residual scheme dependence of the evolution. 

The numerical value for the invariant mass M defined in eq. (|2l|) and calculated 
in the MS scheme with a s (M 2 ) = 0.343 reads 

M = 312 ± 65 cxp ± 7 <f .) ± 14 th MeV . (36) 

This value can be used for comparison of the results of the strange quark mass deter- 
mination obtained from different theoretical calculations and experimental data. If 
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the effective scheme is used for determination of the numerical value for the invariant 
mass M from eq. ( pip the result is slightly different that reflects the residual scheme 
dependence due to truncation of PT series for (3- and 7-functions. 

Results presented in eqs. ( p5| , ^4|) are the main new numerical results of the paper. 



Another new results is a formulation of the problem in an effective scheme and the 
development of all the necessary techniques for phenomenological applications. 

Note that the uncertainty of the final result is smaller than the difference between 
the results of mass extraction from the zero and second moment (first and the last 
lines of Table Q). As we have mentioned before and explain in detail later on we 
do not consider the results obtained from the high moments as reliable which means 
that the theoretical uncertainties of the numerical mass value my' = 111 MeV (from 
the first moment) and m^ = 95 MeV (from the second moment) are expected to 
be much larger than the uncertainty of our optimal choice - the zeroth moment. 
Therefore, the results of all determinations are consistent but the uncertainty of higher 
moments is large. Our view is that after accounting for high dimension condensates 
the second and third determination will change (while the first one does not) and may 
lie within the error bars. Since there is no reason to expect the effect of high-dimension 
condensates to be negligible the requirement that the first and third determination are 
the same implicitly implies that high- dimension condensates are negligible which is 
just an additional assumption (without any justification). Because we do not use these 
moments in our determination we do not quantitatively discuss these uncertainties. 

4.2 Comparison with other results 

Our final result for the numerical coefficient in front of the ml correction in eq. ( [33]) 
agrees with the estimate obtained from the analysis in the MS scheme ||26|| . In the 



present paper we find 20.2 ± 1.8 for the coefficient of the m s correction while the 
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final result of ref. |26| with the conservative estimate of the error bar is 18.1 ± 4.8. 
The final uncertainty of the result obtained in ref. [^J was determined from both 
resummed and FOPT analyses. If only the CIPT analysis in the MS scheme is used 
then one obtains the value 18.1 ± 2.6 which has a smaller uncertainty [p6fl . The 
present value 20.2 ±1.8 results from the analysis in the effective scheme with stricter 
criteria of convergence and more conservative error estimate; for instance, in the final 
result for F^i we have doubled the error for coefficients uo q ^ g from eqs. (|24"|). However, 
the interpretation of the perturbation series in a closed manner (resummation in the 
effective scheme) allows one to reduce the uncertainty of theoretical expressions in 



comparison with the previous analysis in the MS scheme performed in ref. ||26|| . This 
leads to an essential reduction of the theoretical part of the error in the extracted 
numerical value for the strange quark mass. 



The analysis of mass corrections in the MS scheme was also done in ref. p7 | 
where a very accurate account for known PT corrections to the coefficient functions 



of nonPT corrections was given. The direct comparison of our results with ref. [27 



is not simple because different representations have been chosen for the observables. 
Also some approximations have been used which prevents us from directly comparing 
with the results of the present paper at intermediate stages. The final results for the 



mass extracted from the zeroth order moment is rather close to the result of ref. [27 



Still in general the results of refs. |26|. ^7| are rather close to each other concerning 
use of scheme while the change of the scheme was the main reason for our analysis 



as for comparison with ref. [26]. The analysis of the present paper also extends the 
analysis of scale dependence which has been done in ref. [ 27f1 . 

The resulting values of the strange quark mass are close to the earlier estimates 
obtained in refs. ]45], ||| [|7|, H] with the use of less accurate theoretical input formulas 
and less precise experimental data. The recent analysis based on (pseudo) scalar 
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correlators with a thorough parameterization of experimental data gives a value close 



to ours 1 49, 50]. In a lattice calculation of the numerical value for the strange quark 
mass [52] the theoretical input is of nonperturbative nature. The recent results 



obtained in lattice calculations are smaller than our value for m s but still inside the 
error bars as given in eq. (j34|). 

4.3 Choice of moments 

High moments with the weight function (1 — s/M^) k (large values of k) are saturated 
by nonperturbative (infrared, or low-energy) contributions because the perturbative 
region of integration (large energy) is suppressed. On the experimental side this 
means that such moments are saturated by the contributions of low lying resonances. 
In the considered case this is the contribution of the K meson in spin channel and of 
the system i^i(1270) — i^i(1400) in the spin 1 channel. The low-energy contributions 
can be accurately measured that makes the moments obtained from the experimental 
data rather precise. On the theoretical side, within the OPE calculation of the cor- 
relator the infrared sensitivity of high-A; moments means that the PT contribution to 
such moments is suppressed and the moments are saturated by the contributions of 
vacuum condensates of high-dimension operators. Implicitly this is seen in a poor con- 



vergence of the PT series for such moments which is demonstrated in [ FJ| . To obtain 
an accurate numerical value for such moments one has to include vacuum condensates 
of high-dimension operators. However, numerical values for vacuum condensates of 
the local operators with dimension larger than six are completely unknown and are 
usually neglected in the r-system analyses. Therefore, high moments with the weight 
function (1 — s/M^) h have an uncontrollable admixture of high- dimension condensates 
that makes them strongly nonperturbative and, therefore, unreliable for practical PT 
applications. Indeed, the moments (k, 0) obtain contributions from all condensates 
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up to dimension D = 2k + 8 the numerical values of which are unknown for large k. 
Because on the experimental side this contribution corresponds to the contribution 
of low-lying resonance which is not described by PT it seems unjustified to neglect 
the contributions from high-dimension condensates and use only the PT expressions 
(supplemented by condensates only up to dimension 6 which are available) for large k 
moments. In the analysis presented in ref. ]27| the contributions from high- dimension 
condensates (with D > 6) were regarded as an additional theoretical uncertainty. 
Note, however, that the high- dimension condensates contribute quite differently to 
moments (k, 0) and (0, 1) because the integrals corresponding to these moments are 
saturated by contributions of different energy regions. The (0, 1) moments for large 
I are saturated at large energies (of order M T ) and nonPT contributions described 
by the high-dimension condensates within OPE are small. In the leading order only 
one condensate is picked up by integration for a {0,1) moment (assuming that the 
factor (1 — s/M 2 ) 2 is removed from the rate). On the contrary, the (k,0) moments 
for large k are saturated at small energies which correspond to the region of strong 
coupling: these moments are dominated by the resonances. Therefore the (k, 0) mo- 
ments are definitely nonperturbative for large k. Within OPE paradigm this means 
that the total contribution of high- dimension condensates is large compared to per- 
turbative contribution. Indeed, all condensates with dimension up to D = 2k + 8 
contribute to (k, 0) moments and arrange themselves in a way to reproduce contri- 
butions of the low-lying resonances according to the standard phenomenology within 
OPE. This situation is clearly seen in exactly solvable models where all power correc- 
tions are known. In QCD, however, numerical values of high-dimension condensates 
are not known and the total contribution of high- dimension condensates cannot be 
quantitatively analyzed for large k moments though qualitative arguments are rather 
transparent. These reasons forced us to use only the (0, 0) moment as the most re- 
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liable from the theoretical point of view even despite the fact that the experimental 
precision for high k moments is better. 

5 Conclusion 

We have considered the m^-corrections in a QCD-based description of the r-system. 
We use a natural effective scheme well suited for a contour improved perturbation 
theory analysis. The quality of our results is determined by the PT expansions 
of the effective (5- and 7-functions which are the only perturbative objects in our 
analysis. The 7 q -function already shows an asymptotic growth at NNLO while the 
f3- and 7 9 -functions still 'converge' up to this order. In our discussion of the N 3 LO 
terms of PT series for these functions (which depend on the unknown parameter k 3 ) 
we found strong indications for asymptotic growth at this order. This shows that 
the ultimate theoretical limit of FOPT precision is already reached for the set of 
observables in the r-system. This is not yet an actual problem of QCD because of 
the insufficient precision of the experimental data on Cabibbo-suppressed r decays 
especially regarding the differential r decay rate. The experimental situation may, 
however, improve soon. Then our procedure of using the effective scheme description 
of the r-system opens the possibility of high precision tests of QCD independently 
of an explicit convergence of PT series in the MS scheme. In this field an effective 
scheme approach can show its real power because then the main source of theoretical 
uncertainty, the relation of the internal mass parameters m q , m g to the MS parameter 
m s , will be eliminated. Note that for a QCD test within our approach it is necessary 
to relate four r-observables to each other because the three parameters a, m g and m q 
have to be fixed. 

Our result for the strange quark mass is m s (M^) = 130 ± 27 cxp ± 9 th MeV where 
we have combined the pure theoretical error which is basically determined by the 
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truncation of PT series for effective RG functions and the error due to the strange 
quark condensate into one number linearly. After running to the scale 1 GeV we 
obtain m s (l GeV) = 176 ± 37 exp ± 13th MeV. This is consistent with the previous 



results where the resummation was done in the MS scheme [26|. The large part of 
the theoretical uncertainty of our result for the PT coefficient of m 2 s term comes from 
re-expressing the effective quantities in terms of the MS scheme parameters. The 
advantage of the procedure presented here is that the estimate of the accuracy is 
not based on the decomposition of the result into terms coming from corrections to 
the correlator in the MS scheme. This decomposition seems unnatural in resummed 
perturbation theory as it introduces an additional uncertainty related to the conver- 
gence of the series for the correlator in the MS scheme every term of which is given 
by a closed expression resulting from the resummation along the contour. Within the 
effective scheme approach all sources of uncertainty are collected into the effective 
(3- and 7-functions which are the only PT quantities entering into the analysis that 
provides a solid ground for estimating the accuracy of theoretical expressions. 
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Appendix 

We give the expressions for the effective (3- and 7-function coefficients in terms of the 
standard MS scheme results. In this Appendix the coefficients of (3- and 7-functions 
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without upper index stand for MS quantities while the coefficients with upper index 
"eff" denote the coefficients of the effective functions. 

The coefficients of the effective /^-function (|8]) are given by 

Ptf = A), P'f = Pi, 

/3f = (3 2 -k 1 /3 1 + (k 2 -k 2 1 )(3o, 

(3f = p 3 -2^ 2 + ^+ (2k 3 -6fc 2 fci+4fc?)A>- (37) 

For the effective 7-function coefficients of eqs. fl9j|T0l) in terms of the MS scheme (3- 
and 7-function coefficients we find 

1 



eff _ _ eff 

= 71 _ K H0 1" 2 re n( 

!.„/., , 1 



7n0 = 70, 7nl = 71 - hlQ + -KoPo 



7n2 = 72 - 2fei7i + (~k 2 + 2fci)7o + -fcnoA + (-&l&n0 + &nl - ^no) A), 

7^3 = 73 - 3fc l72 + (-2A: 2 + 5fc?) 7 i + (_fc 3 + 5^*! - 5fc?) 7o 

+ 2 k n0@2 + f-^l&nQ + ~ \ k loj Pi ( 38 ) 

+ (^—k 2 k n Q + -klk n0 — Skxkni + -Aa&^ + T^™ 2 — T^™ 1 ^™ 2^™°) 

Here fc n j with n = q,g stand for the coefficients of eq. @ if n = g (k gQ = 5/3) and 
for the coefficients of eq. (0) if n — q (k q0 = 7/3). Numerical values for all necessary 



coefficients in the MS scheme have been collected in |12| where further references to 
original papers can be found. 
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